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ABSTRACT 

Aims. The main goal of our project is to investigate the spatial distribution of different stellar populations in the Magellanic Clouds. 
The results from modelling the Magellanic Clouds can also be useful for simulations during the Gaia mission preparation. 
Methods. Isodensity contour maps have been used in order to trace the morphology of the different stellar populations and estimate 
the size of these structures. Moreover, star density maps are constructed through star counts and projected radial density profiles 
obtained. Fitting exponential disk and King law curves to the spatial distribution allows us to derive the structural parameters that 
describe these profiles. 

Results. The morphological structure and spatial distributions of various stellar components in the Magellanic Clouds (young and 
intermediate age stars, carbon stars), along with the overall spatial distribution in both Clouds, are provided. 

Key words, galaxies: Magellanic Clouds 



1. Introduction 

The structure of the Mag ellanic Clouds has bee n a su bject 
of study for many years. Ivan der Marel & Cionil (l200ll) and 
Ivan der M arel ( 200 1 ) determined the viewing angles of the LMC 
and constructed a near-IR star count map, demonstrating that the 
LMC is intrins i cally elongated. Many articles of the past decades 
dBruckl 119801: iDottori et al.1 Il996t iGardiner & Hatzidimitrioul 
also discuss the different area distributions of the vari- 
ous stellar populatio ns in the MCs. iMaragoudaki et all (1 1998b . 
\ IZaritskv et"aiT(l2000l) . and IMaragoudaki et all d2001l) show that, 
in both LMC and SMC, the spatial distribution of the stellar 
content tends to become more regular and ordered as its age 
increases. Young stars are found to be results of bursts of star 
formation, although there is some disagreement on the reasons 
for the triggered star formation. 

Consistent with the previous works are the results from pa- 
pers based on infrared surveys. Counts of objects towards the 
Magellanic Clo uds from th e DEN IS near-infrared survey have 
been studied bv lCioni etal . (2000), who used colour-magnitude 
diagrams (CMD) to distinguish between three groups of objects 
with different mean ages. The spatial distribution of the three age 
groups is found to be quite different: the youngest stars exhibit 
an irregular structure, while t he older stars are smoot h ly and reg- 
ularly distributed. Mor eover, iNikolaev" & Weinberg! (|2000) and 
Gonidakis et ak| (120091) have investigated the spatial distribution 
of the LMC and SMC stellar components, respectively, from 
2MASS data. They present morphological analysis of the MCs 
colour-magnitude diagram, identifying the different populations 
and estimating the ir projected sp atial distributions of various 
stellar populations. Bica et all (120081) propose a scenario where 
the present-day, angular-averaged, large-scale structures of both 
Clouds appear to behave as tidally truncated systems (which is 
expected, since they are believed to be Milky Way satellites), 



characterized by well-defined core and halo substructures, even 
if the Clouds are not spherical systems. Thus, they have under- 
gone severe tidal perturbation when the last dynamical and hy- 
drodynamical i nteraction between the Clouds took place (about 
200 Myr ago. lBekki &"Chibal (12007b ). The older LMC and SMC 
star clusters, on the other hand, appear to be distributed as an 
exponential disk. This distribution is possibly reminiscent of the 
Clouds' structure prior to the last interaction. 

Hopefully, some of the questions regarding the Magellanic 
system can be answered by ESA's cornerstone mission Gaia. It 
is primarily an astrometric mission, whose goal is to create the 
largest and most precise three-dimensional map of our Galaxy 
by providing unprecedented positional and radial velocity mea- 
surements for about o ne billion stars in our G alaxy and through- 
out the Local Group dPerrvman et al.ll200ll) . The detection of 
stars is expected to be complete to V=20 mag (astrometry + 
photometry), but the radial velocities will not be measured for 
stars fainter than V=17.5 mag. The astrometric precision should 
be 25 //as at V=15 mag and the distances will be known with a 
precision better than 10% for hundreds of millions of stars. Due 
to the spatial resolution nearby galaxies will be resolved in stars. 
Thus Gaia data will improve estimation of rotational parallaxes 
for Local Group galaxies, kinematic separation of stellar popu- 
lations, galaxy orbits, and cosmological history, all very interest- 
ing topics of our nearest galaxies. The galaxies that will benefit 
most from the Gaia survey will be the Magellanic Clouds, since 
they are nearest to our Galaxy, where millions of stars are ex- 
pected to be detected by the Gaia instruments. 

The main goal of our project is to obtain the spatial distri- 
bution of several stellar components in these galaxies, combin- 
ing the available data sets of various stellar populations. As the 
nearest neighbours of our Galaxy, the Magellanic Clouds are the 
most important targets with a large number of stars observed by 
Gaia. The results of this investigation will also be used to im- 
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Table 1. Age groups 





Age 


magnitude 


colour 


A 


<0.1 Gyr 


11 <U< 16 


-1.5<U-V<-0.6 


B 


0.1-0.3 Gyr 


16<U< 18 


-1.3<U-V<0.2 


C 


0.3 - 0.9 Gyr 


18 <U<21 


-0.8<U-V< 1.2 


D 


> 1 Gyr 


17 < V < 20 


0.4 < V-I < 1.8 


E 


<0.1 Gyr 


K< 15 


-0.5 < J-K < 0.2 


F 


> 1 Gyr 


13.5<K<15.5 


0.2 < J-K< 1.2 



prove the Gaia Universe Model (Robin et al. 20091), a set of al- 
gorithms used by the data generators of the Gaia simulators to 
generate simulated data in the framework of Gaia data reduction 
preparation. Our contribution to the Universe Model consists in 
implementing a model for the galaxies, resolved in stars, as seen 
by Gaia. In this paper we present our first results from mod- 
elling the density distribution of the Magellanic Clouds based 
on various available data sets that are appropriate as described in 
section 2. In section 3 the methodology we follow is explained, 
and in section 4 the results of this work are presented. Finally 
our concluding remarks are gathered in section 5. 



2. The data 

To obtain the spatial distribution we used archive data. We se- 
lected all-sky surveys and dedicated catalogues, which are ho- 
mogeneous, have a good coverage of the galaxies, and are deep 
enough. We are also interested in the distribution of more spe- 
cific objects, such as carbon stars, or of extended objects like star 
clusters and associations. 



2.1. The Magellanic Clouds Photometric Survey (MCPS) 

The Magell anic Clouds Photometric Survey (Zaritsk y et al.l 
I20021 l200l offers very rich and complete catalogues of the 
Magellanic Clouds, but the area they cover is s omew hat lim- 
ited. The M CPS for the SMC by IZaritskv et all d2002l) and for 
the LMC bv lZaritskvetal.1 d2004l) are catalogues of the U, B, V, 
and I stellar photometry of only the central 18 deg 2 area of the 
SMC and of the central 64 deg 2 area of the LMC. The incom- 
pleteness becomes significant at a magnitude fainter than V < 20. 
This is the reason we are limiting our study to stars brighter than 
V=20fl The data are combined with 2MASS and Deep Near- 
Infrared Southern Sky Survey (DENIS) catalogues to provide, 
when available, U through K, data for the stars. 

2.2. The Two-Micron All Sky Survey (2MASS) 

2MASS has uniformly scanned the entire sky in three 
near-infrared bands to detect and characterize point sources 
brighter than about 1 mJy in each band, with a signal-to- 
noise ratio greater than 10. With 2MASS we can choose 
as large area around the Magellanic Clouds as is necessary. 
However, the limiting magnitude is somewhat bright with 
J<15.8, H<15.1, and Ks<14.3 mag dSkrutskie et alj |2006). The 
data used here are from the 2MASS All-Sky Point Source 
Catalog, at IPAC Infrared Science Archive (IRSA), Caltech/JPL 
(http://irsa.ipac.caltech.edu/applications/Gator). 



We used the c riteria for the main sequence and red giant stars 
bv lGavras[ (2003) (see TableQ}. These also agree with the criteria 
used bv lNikolaev & Weinberg] d2000h . 

2.3. Carbon stars 

We used catalogues of carbon stars in the LMC and SMC from 
objective-prism plates, taken with the UK Schmidt Te l escope . 
The LMC catalogue is produced by Kontizas et al. (H 001), 
and the SMC catalogues are fr om iRebeirot et al. (1 19931) and 
iMorgan & Hatz idimitrioul d 19951) . 

Hess diagrams from the various data sets are produced. A 
Hess diagram shows the relative density of occurrence of stars 
at different colour-magnitude positions in the colour-magnitude 
diagram. The colour-magnitude range from the data sets we 
use was divided in a grid with different numbers of cells for 
each data set. The size of the cells should provide a fine 
grid and at the same time contain enough stars for the statis- 
tics. The diagrams are shown in FigQ] We divided the stel- 
lar content of both Magellanic Clouds into several age groups, 
matching features of the CMD with i sochrones ob t ained from 
http://stev.oapd.inaf .it/cmd based on Mar igo et al.l d2008l) and 
iBertelli et all (ll994l) .The isochrones for 10 Myr, 100 Myr, 300 
Myr, and 1 Gyr are displayed on top of the contours in the Hess 
diagrams. The selection criteria we use are summarized in Table 



1 This is convenient since stars fainter than 20th magnitude will not 
be detected by Gaia because of its observational limit. 



3. Methodology 

To determine approximately the shape and the distribution of the 
galaxies we used two approaches. 

1 . We performed star counts with a rectangular grid. The size of 
the grid cell was different for the various stellar populations. It 
was chosen in such a way as to include a sufficient number of 
stars and also to provide enough detail for the isodensity contour 
maps. The contours in these maps trace areas with equal stellar 
density. They show the overall shape of the galaxies and offer an 
initial insight into their spatial distribution. 

2. We also used radial density profiles (RDPs) to obtain the ac- 
tual distribution of the galaxies. The RDPs correspond to the 
projected radial number density of objects contained in concen- 
tric rings around the LMC and SMC centroids. The underlying 
assumption for this kind of analysis is that the structures should 
present an important degree of radial symmetry. This is not al- 
ways the case for the Magellanic Clouds, but RDPs can still be 
used as probes of the radial distribution of the various objects 
averaged over all azimuthal directions, hence of the large-scale 
structure. We divided the area around each galaxy in concentric 
annuli of increasing steps outwards and we counted the number 
of stars within each annulus. The number of stars per unit area 
in each annulus gives us the stellar density, corresponding to the 
mean distance from the centre of the galaxy. The thus derived 
density profiles for each age and/or object group are compared 
with existing theoretical models 

For the centres of the Magellanic Clouds, we adopt the 
following coordinates given by the SIMBAD Astronomical 
Database (http://simbad.u-strasbg.fr/simbad) - for LMC RA: 
05 A 23'", Dec: -69" 45 m and for the SMC RA: 00' 1 52 m , Dec: 
-12 d 48'" 

A conventional first approach to describing the spatial dis- 
tribution of various stellar populations in the Magellanic Clouds 
from the RDPs is to use an exponential-disk profile, 



fir) = f 0D x e 
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M. K. Belcheva et al.: Spatial distribution of stellar populations in the Magellanic Clouds 



3 






-0.5 0.0 0.5 1.0 1.5 2.0 
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Fig. 1. Colour-magnitude diagrams of the LMC from the MCPS (left and middle panel) and 2MASS (right panel). The density levels 
are logarithmic. Isochrones for 10 Myr, 100 Myr, 300 Myr and 1 Gyr are displayed on top of the contours. 



where ho and foo are the scale length and the central density of 
objects, respectively, and r is the distance from the centre. As a 
secon d approac h we use a King-like profile based on the King 
law {King 1962), which is often used to describe the distribution 
of globular clusters, but also applies to dwarf spheroidal galax- 
ies, 



f(r) = fa 



where r, and r c are the core and tidal radii, respectively, and 
fox the central density of objects and r is the distance from the 
centre. 

The best-fitting profiles are found by performing the 
Levenberg-Marqwardt least-squares fit to the considered func- 
tions. This fit was performed with IDL (Interactive Data 
Language, http://www.ittvis.com/) with the use of th e specially 
developed procedure MPFITFUN (lMarkwardlll2009l) . 

4. Results 

Figure [21 |4j [6j and[8]show the isodensity contour maps obtained 
from the MCPS, 2MASS and carbon star catalogues. The spa- 
tial distribution of the various age groups is as expected. The 
older stars have a more regular and smoother appearance, while 
younger stars form fragmented and less symmetric structures 
(Cioni et al. l2000UZaritskv et al.l2000tlMar agoudak i et al.ll998l 
l200lHGonidakis et alJl2009l) . 

The RDPs obtained from MCPS, 2MASS, and the carbon 
stars are presented in Figure [3] [5] |7l and [9] The best fit of both 
models is displayed in each figure. In all figures, except Figure 
|9l error bars are not shown, because the Poisson errors are com- 
parable in size to or smaller than the symbols. 

The structural parameters of the LMC and SMC obtained 
by fitting RDPs from the adopted data sets in the various cata- 
logues are listed in Table|2] Columns (2) and (3) contain the cen- 
tral density of objects foo and the scale length ho respectively. 
Columns (4)-(7) contain the central density of objects fox, the 
core radius r c , the tidal radius r, where available, and the con- 
centration parameter c p . The parametric value of the tidal radius 
r, obtained for stars from 2MASS and the MCPS was not reli- 
able enough, because of the relative incompleteness of 2MASS 
and the small area that is available from the MCPS. The carbon 



stars, on the other hand, represent a sufficiently complete sample 
of this particular group of objects as numbers, spatial distribu- 
tion, and magnitude range. Therefore we can use these objects 
as prototypes of the major dynamical system in both Clouds. 

Figures [3^ and [5^ display a very clumpy distribution for 
the brightest stellar population. Considering that the stars in 
this age group are in concent rations, i.e. stellar complexes 
jMaragoudaki et alJI 19981 1200 lb . then the very wide spread is 
expected and explained. 

The RDPs for the SMC from MCPS show that both the King 
and the exponential-disk models can be used to describe the stel- 
lar distribution. However, as can be seen from Fig. [5] if the dis- 
tribution of the youngest stars is almost equally well fitted by 
both, the oldest stars prefer the King law. Such a difference is not 
observed in the LMC, where the exponential-disk model seems 
slightly better. This is possibly due to the small area available, 
which is dominated by the younger stars. A similar behaviour is 
observed in the RDPs from 2MASS. Both the younger and the 
older stars in the LMC are distributed on an exponential disk, 
while the radial density distribution of the SMC is better de- 
scribed by the King profile. 

Additionally, the residuals below each profile reveal the pre- 
ferred model and show the dominant model for each part of the 
galaxy. 

The spatial distribution of stars from 2MASS and of the car- 
bon stars in the LMC (see Figs.[6]and[8]) suggests that there are 
two subsystems, whose major axes are almost perpendicular to 
one another. The same behaviour has recently been detected by 
iBica et alJ(l2008l) for the star clusters in the LMC. T his effect has 
previously been des cribed bv lKontizas et al.l (Il990h and later by 
Dotto ri et alj (1996), who found that the young clusters, mainly 
occupying the bar, are rotated with respect to the older clusters. 

5. Discussion & conclusions 

The structural parameters of the LMC and SMC from the 
adopted data sets in the various catalogues (MCPS, 2MASS, and 
carbon stars) are listed in Table [2] Neighter the exponential-disk 
nor King models are always very explicitly determined as the 
dominant dynamical model. However, we can adopt a criterion 
for the acceptance of one or the other model. The values given 
by the x 2 test allow us to assume one of the models as more 
appropriate for describing the distribution of the data. A second 
criterion could be the age of the stellar population, assuming that 
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Fig. 2. Isodensity contour maps of the LMC stars from the MCPS with various ages. Top row - age groups A and B; bottom row - 
age groups C and D. 

Table 2. Structural parameters measured from the RDPs with the exponential-disk and King profiles 



RDP 


fi)D 


h D 


foK 


r c 


r t 


c p 




[Obj. deg- 2 ] 


[deg] 


[Obj. deg- 2 ] 


[deg] 


[deg] 




MCPS 














LMC A 


(1.4 ± 0.9) x 10 3 


2.6 + 0.3 


(1.2 + 0.6) x 10 3 


2.1+0.2 






LMC B 


(18+ l)x 10 3 


1.33 + 0.07 


(16+ l)x 10 3 


1.1+0.06 






LMC C 


(96 ± 4) x 10 3 


1.44 + 0.06 


(83 ± 4) x 10 3 


1.2 + 0.07 






LMC D 


(260 + 3) x 10 3 


1.52 + 0.02 


(225 ± 4) x 10 3 


1.3+0.03 






SMC A 


(7 ± 0.6) x 10 J 


0.57 ± 0.04 


(6.7 ± 0.5) x 10 3 


0.4 ± 0.03 






SMCB 


(64 ± 3) x 10 3 


0.48 ± 0.02 


(62 ± 2) x 10 3 


0.4 + 0.01 






SMCC 


(183 ±4) x 10 3 


0.64 ± 0.01 


(177 ± 6) x 10 3 


0.6 ± 0.02 






SMCD 


(259 ± 11) x 10 3 


1.13 + 0.06 


(918 ± 247) x 10 3 


1.9 + 0.2 






2MASS 














LMC E 


(7.7 ± 0.2) x 10 2 


2.01+0.05 


(11.1 ± 1.0) x 10 2 


2.89 ± 0.2 






LMC F 


(41.3+ 0.5) x 10 3 


1.89 + 0.02 


(36.3 ± 0.4) x 10 3 


1.58 + 0.02 






SMCE 


(22.6 ± 1.0) x 10 2 


0.60 ± 0.01 


(23.8 ± 0.7) x 10 2 


0.40 ± 0.01 






SMCF 


(28.7 ± 0.7) x 10 3 


1.03+0.02 


(26.3 ± 0.3) x 10 3 


0.82 + 0.01 






Carbon stars 














LMC 


368 ± 22 


1.69 + 0.07 


514+17 


3.3 + 0.1 


9.7 + 0.2 


0.46 


LMC bright 


72 ±6 


0.98 ± 0.05 


71+7 


0.9 + 0.1 


12 + 3 


1.12 


SMC 


1050 ± 40 


0.69 ± 0.02 


1100 + 39 


0.85 ± 0.04 


5.2 ± 0.3 


0.78 



the very young population might be distributed on an exponen- 
tial disk-like structure, even if it does not fit the whole galaxy. 

It was shown that there is a system of young objects more 
concentrated in the central region of the LMC with a position 
angle almost perpendicular to the other system. Of course the 
inclination of the LMC is small (almost face-on), allowing the 
difference to be seen. On the other hand, the SMC has a high 
inclination to offer clear evidence of such a characteristic. The 
two systems in the LMC have been previously mentioned, and 



are explained as one result of the interaction between the Milky 
Way galaxy and the Magellanic Clouds. 

From the carbon stars in the LMC, it was shown that two 
different systems exist with a core radius of 3.3 + 0.1 deg for the 
faint and 0.9+0. 1 deg for the bright carbon stars. An exponential 
model also shows such a difference in the scale height, 1 .69 + 
0.07 deg and 0.98 + 0.05 deg, respectively, revealing a smaller 
central system of more massive carbon stars. The young clusters 
are also found in a smaller central system, than for the old ones 
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dBica et alJl2008h . We did not observe such a segregation for the 
SMC carbon stars. 

Although the Magellanic Clouds are assumed to be irregu- 
lar galaxies, it also seems that the older populations appear to 
behave as tidally truncated systems, even if they do not show 
obvious radial symmetry. From Table [2] we can see that the pa- 
rameters of the fitting for the young stellar populations are in- 
conclusive. However, this is no surprise, because even if the bulk 
of the stars show some radial symmetry, this is not the case for 
the young bright stars. Both their contour maps and their radial 
density profiles do not support such an assumption. This is an- 
ticipated since clumping of star-forming regions dominates their 
distribution. 

The results of this investigation were used to provide data 
and model parameters for Gaia simulations and in particular for 
the Gaia Universe Model. 
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Fig. 3. RDPs for the subsets of LMC stars from the MCPS, fitted with exponential-disk (solid line) and King profiles (dashed 
line). Top row - age groups A and B; bottom row - age groups C and D. Error bars are not shown, because the Poisson errors are 
comparable in size to or smaller than the symbols. Below each profile are the residuals from the fit. 
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Fig. 5. RDPs for the subsets of SMC stars from the MCPS, fitted with exponential-disk (solid line) and King profiles (dashed 
line). Top row - age groups A and B; bottom row - age groups C and D. Error bars are not shown, because the Poisson errors are 
comparable in size to or smaller than the symbols. Below each profile are the residuals from the fit. 




Fig. 6. Isodensity contour maps of the Magellanic Clouds from 2MASS. Top row - age group E; bottom row - age group F. Left - 
LMC, right - SMC. 
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Fig. 7. RDPs for the subsets of stars from 2MASS, fitted with exponential-disk (solid line) and King profiles (dashed line). Top 
row - age group E; bottom row - age group F. Left - LMC; right - SMC. Error bars are not shown, because the poisson errors are 
comparable in size to or smaller than the symbols. Below each profile are the residuals from the fit. 



M. K. Belcheva et al. 



: Spatial distribution of stellar populations in the Magellanic Clouds 



11 




Fig. 9. RDPs for the carbon stars, fitted with exponential-disk (solid line) and King profiles (dashed line). Left - LMC; middle - 
LMC bright; right - SMC. Below each profile are the residuals from the fit. 



